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Fig. 4. (A) Mean number of mobbing alarm calls produced by 15 fairy- 
wren groups in response to different model types: cuckoo (CK), currawong 
(CR), honeyeater (HE), snake (SN), and sparrowhawk (SP). (B) Mean time 
spent mobbing the cuckoo model (<0.5 m from the model) by small (n = 27 
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groups) versus large (a? = 5 groups) groups. (C) Mean number of whining calls 
produced in response to the five model types. (D) Mean number of individuals that 
approached playbacks {n = 20 each) of fairy-wren mobbing alarm calls, a control 
sound, and fairy-wren whining alarm calls. Error bars denote the standard error. 



control sound (a parrot call, GLMM: x\ = 68.05, 
P < 0.0001 ; Fig. 4D). Once mobilized, the strength 
of these defenses increases with group size. Large 
groups were more vigilant around their nest 
(GLMM: x^i = 8.03, P < 0.004), spent more time 
mobbing the cuckoo than smaller groups (Kruskal 
Wallis test: x^i = 5.42, P = 0.02; Fig. 4B), and 
ultimately were less likely to be parasitized. Thus, 
superb fairy-wrens possess cuckoo-specific nest 
defenses, which are enhanced by helper con- 
tributions and which can explain the lower 
parasitism rates experienced by large groups. 

Our findings show a pronounced association 
between avian brood parasitism and cooperative 
breeding in birds, on a global scale. Our field data 
suggest that a two-way process underpins this 
relationship. On the one hand, brood parasites 
can gain a fitness advantage by preferentially ex- 
ploiting the superior care provided by coopera- 
tively breeding groups. On the other hand, the 
genetic relatives of offspring raised by coopera- 
tively breeding families potentially gain fitness 
from the superior defenses that the extended fam- 
ily collectively mounts against brood parasites. 
Defense against brood parasitism is therefore an 
important kin-selected fitness advantage associ- 
ated with cooperative breeding [see also (7(5)]. In 
superb fairy-wrens, we have shown this two-way 
process at work, but in other cooperatively breed- 
ing hosts, especially those with a less protracted 
coevolutionary relationship with brood parasites, 
only the first part of the process may be evident. 
The challenge remaining for fixture work is to 
determine the extent to which brood parasites 
have influenced the biology of cooperatively breed- 
ing species. 
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and C57BL/6N became a second major source of 
C57BL/6 mice. Laige-scale projects use different 
C57BL/6 substrains, including the Intemational 
Knockout Mouse Consortium (IKMC), which uses 
C57BL/6N embryonic stem (ES) ceUs (2, 5), and 
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the Mouse Genome Sequencing Consortium and 
Allen Brain Atlas, which use C57BL/6J (4, 5). 
Although most behavioral measurements to date 
have been carried out in the C57BL/6J substrain, 
the IKMC will shift emphasis to the C57BL/6N 
substrain. Although behavioral differences have 
been noted in C57BL/6 substrains, their genetic 
basis has not been elucidated (6, 7). Therefore, it 
becomes important to understand the genetic causes 
of phenotypic differences between C57BL/6N and 
C57BL/6J. 

Drug addiction is characterized by loss of con- 
trol over drug consumption and behaviors asso- 
ciated with drug seeking. Even though drugs of 



REPORTS 

abuse fall into different pharmacological cat- 
egories, they all act on the mesolimbic reward 
pathway in the brain (5). Long-term structural 
and fiinctional changes in neuronal circuitry are 
thought to be a key feature of addiction (P). We 
investigated in detail a difference in cocaine re- 
sponse between the C57BL/6N and the C57BL/6J 
mouse substrains. 

We characterized C57BL/6N and C57BL/6J 
animals for locomotor response to cocaine. Be- 
cause of the nonlinear nature of drug response, 
we carried out testing at multiple doses. C57BL/6N 
had a 45% (1 SD) lower acute response to cocaine 
and methamphetamine at multiple doses (Fig. 1, A 
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Fig. 1. Acute and sensitized psychostimulant response, measured by 
locomotor hyperactivity, is lower in C57BL/6N (B6N) substrain than in 
C57BL/6] (86]). (A) Baseline locomotor velocity data of B6] (blue) and B6N 
(green) mice were measured for 30 min, and then mice were injected with 
cocaine (20 mg/kg, red arrow) and recorded for 60 min further. (B) Cocaine 
dose response for B6] and B6N. (C) Decreased methamphetamine response in 
B6N. Animals were injected with methamphetamine (2 mg/kg, red arrow) and 
recorded for 3 hours. (D) Methamphetamine dose response for B6] and B6N. 
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Mice were treated with 1, 2, and 4 mg/kg of the drug. Average response 60 min 
post injection is shown. (E and F) B6N have a lower sensitized response to cocaine 
at 10 mg/kg (E) but respond similarly at 15 mg/kg dose (F). Mice were measured for 
baseline and postinjection response for 1 hour and injected with saline (S) or 
cocaine (C). The 30-min postinjection response is shown. All data are shown as 
mean ± SEM; the number of animals in each group is indicated. *P < 0.05,**P < 
0.01,***P < 0.001,****P < 0.0001 from Tukey post hoc test following two-way 
analysis of variance (ANOVA) [(B) and (D)] or pairwise t test [(E) and (F)]. 
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to D, and fig. SI). In psychomotor sensitization, 
repeated administration of the stimulus elicits a 
progressively larger behavioral response. It is reg- 
ulated by experience-dependent neuronal plasticity 
and thought to be a key event leading to addiction 
{10, 11). At a 10-mg/kg dose, C57BL/6J sensitized to 
cocaine much more efficiently than C57BL/6N 
(conparc day 5 to 10, Fig. IE), although both strains 
sensitized to similar levels at a higher 15-mg/kg 
dose (Fig. IF). 

To determine which genetic loci contribute 
to reduced cocaine response seen in C57BL/6N, 
we carried out quantitative trait locus (QTL) analysis 
(fig. S2). Parental strains, Fi, and F2 animals 
were tested concurrently. C57BL/6N had an 
~1-SD lower response to cocaine. The Fi animals 
exhibited a cocaine response similar to that of 
C57BL/6N, indicating dominance of the C57BL/6N 
allele, and the F2 progeny were intermediate in 
response when compared with the parental 
strains (Fig. 2A). These differences were seen 
in all of our measures and in both sexes (figs. S3 
and S4). 

To perform QTL analysis, we identified poly- 
morphic markers between C57BL/6J and C57BL/ 
6N (supplementary materials, fig. S5, and tables 
SI and S2). QTL analysis for cocaine response 



yielded a single QTL on chromosome 11 [log of 
odds (LOD) = 6.8] (Fig. 2C). This peak is highly 
significant on the basis of permutation tests (Fig. 
2C), is specific to cocaine response, and was not 
seen for baseline activity (fig. S6). The peak link- 
age occurs at marker rs 1348 1014 for 30 min, 
60 min, and net response (Fig. 2D), and the geno- 
type effect plot indicates that the B6N/B6N allele 
elicits a lower response than the B6J/B6J allele. 
The Fi response is similar to that of C57BL/6N, 
implying dominance of B6N allele (Fig. 2E). This 
QTL accounts for 11% of the total phenotypic 
variance (genetic, environmental, and error) and 
61% of the genetic variance (fig. S7). The QTL 
support interval translates to a 22-Mb interval 
between 35 to 57 Mb of chromosome 11 (72) (sup- 
plementary materials and fig. S8). 

Phenotypic differences between related strains 
may occur because of three possibilities: residual 
heterozygosity at the time of separation, genetic 
contamination after separation, or genetic drift. 
Genotyping indicated that inbreeding was essen- 
tially complete and that there was no genetic con- 
tamination or residual heterozygosity in C57BL/6N 
and C57BL/6J substrains {13). Because C57BL/6N 
and C57BL/6J diverged 62 years ago, there may 
be a limited number of polymorphisms accumu- 



lated through genetic drift that accounts for the 
phenotypic difference. To identify the causative 
mutation, we performed whole-genome sequenc- 
ing mns and combined the data with published 
and unpublished C57BL/6NJ sequences for 99.8- 
fold coverage of chromosome 11 (table S3) {14). 
We reanalyzed the data sets with identical param- 
eters (fig. S9 and tables S4 and S5) and found a 
single nonsynonymous single-nucleotide polymor- 
phism (SNP) in the QTL interval (Fig. 3A and fig. 
SIO). No indels or stmctural variants cause protein- 
coding changes (Fig. 3, B and C). This nonsyn- 
onymous SNP changes G to A at 46,036,117 base 
pair (bp) of chromosome 1 1 (mm9) and produces 
a serine-to-phenylalanine missense mutation at po- 
sition 968 (S968F) in Cyfip2 (Fig. 3D). This variant 
had a high Polyphen2 score of 0.957/1 and was 
predicted to be damaging (75) (fig. Sll). CYFIPl 
and CYFIP2 are highly conserved proteins im- 
plicated in Fragile X-mental retardation protein 
(FMRP) function and regulation of actin poly- 
merization through the WAVE regulatory complex 
(WRC) (7(5, 17). Cyfip2 is widely expressed through- 
out the brain, whereas Cyflpl is more restricted 
(fig. SI 2). FMRP is the most common cause of 
mental retardation in humans, and WAVE com- 
plex members have been shown to regulate neu- 
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Fig. 2. QTL on chromosome 11 regulates cocaine response. (A) Hyper- 
activity after intraperitoneal injection of cocaine (20 mg/l<g) was quan- 
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single highly significant QTL on chromosome 11. The significance thresh- 
olds are established through 100,000 permutation tests. (D) The chro- 
mosome 11 QTL is significant for multiple cocaine response measures. 
(E) Genotype effect plot at marl<er rsl3481014 shows the B6N allele 
has lower response than B6] allele, and the FI is intermediate. Data are 
mean ± SEAA. 
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ronal connectivity and behavior in mice and 
Drosophila (18-23). 

Because many C57BL/6 substrains are read- 
ily available with information about the time of 
separation from the founder colony, we con- 
structed a phylogenetic timeline of the S969F 
variation (Fig. 3, E and F). The Cyfip2 polymor- 
phism was fixed in the C57BL/6N colony some- 
time between 1961 and 1974 and is present in 
most commercially available sources of C57BL/6N, 
including Charles River (NCrl) and laconic (NTac) 
(Fig.3F). 

The Cyflp gene family is highly conserved in 
metazoans, and a multiple-sequence alignment 
indicated that CYFIP2 S968 is 100% conserved 
in orthologs from Caenorhabditis elegans to hu- 
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mans (fig. S13A). The crystal structure of CYFIPl 
has been solved as part of a pentameric 400-kD 
WRC (16). Molecular modeling of the S969F mu- 
tation using the WRC homology model revealed 
that the mutation leads to substantial steric hindrance 
with adjacent residues (figs. S13B and S14). We 
measured protein stability after cycloheximide 
treatment to estimate the relative half-life of wild- 
type and mutant CYFIP2. The wild-type and mutant 
proteins had half-lives of 8.5 and 2.8 hours, re- 
spectively, indicating that the S968F mutation 
greatly destabilizes CYFIP2 (Fig. 4A) (supple- 
mentary text). Stably transfected HEK293 cell 
lines expressing equivalent levels of the wild- 
type and mutant CYFIP2 proteins were used 
for biochemical and proteomic analyses. Quan- 
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titative co-immunoprecipitation (IP) assays showed 
no difference in interactions between wild-type 
and mutant CYFIP2 with NAPl , WAVEl , ABI2, 
and HSPC300, the core components of the WRC 
(Fig. 4B). We did not observe differences in our 
proteomic analysis (fig. SI 5) and saw no interaction 
of wild-type or mutant CYFIP2 with FMRP, 
FXRl, or FXR2 under our conditions. 

To generate a second mutant allele of Cyfip2, 
we produced mice using ES cells carrying the 
"knockout firsf (null) allele from the IKMC (3) 
and maintained them in C57BL/6N background 
(fig. SI 6). The homozygous knockouts die within 
a few hours of birth; however, the heterozygous 
mice are viable with no overt phenotypes. Because 
the ES cells used in IKMC are of C57BL/6N ori- 
gin, the wild-type mice carry the mutant C57BL/6N 
aUeles of Cyfip2 (designated Cyfipf^^^^\ and the 
heterozygotes carry one copy of the C57BL/6N 
allele (designated Cyfip2^^^^~). We compared 
mice heterozygous for the knockout allele with 
the homozygous WT mice for cocaine response 
(Fig. 4C). As expected the C57BL/6N mice have 
lower cocaine response than C57BL/6J mice at 
all doses (Fig. 4C). The Cj^^^^^^^^^ mice also 
have low response to cocaine, similar to that of 
C57BL/6N; however, heterozygous Cyfip2^^^^- 
mice have higher cocaine response than homo- 
zygous Cyfip2^^^^^^^ mice at all three doses 
(Fig. 4C). The deletion of one mutant allele 
alleviates the severity of the phenotype seen with 
two mutant alleles and explains the dominant 
phenotype of the Fl progeny and the genotype 
effect plot of linked QTL maker seen previously 
(Fig. 2, A and E). The sensitized response is 
also higher in heterozygous mice (Cyfip2^^^^~\ 
confirming that CYFIP2 regulates acute and sen- 
sitized response to cocaine (Fig. 4D). We did not 
see any changes in cocaine metabolism between 
the two strains (fig. SI 7). 

Next, we explored the fiinctional effect of the 
CYFIP2 S968F mutation. Drug-induced struc- 
tural plasticity is thought to be key in addiction. 
We measured and classified dendritic spine mor- 
phology of medium spiny neurons in the nucleus 
accumbens. C57BL/6N have a lower number 
of total spines with significantly lower spines of 
each class (Fig. 4E). Because dendritic spines 
are the major sites of excitatory postsynaptic cur- 
rent, we predicted a deficit in excitatory gluta- 
matergic signaling in C57BL/6N. We measured 
the frequency and the amplitude of mini-excitatory 
postsynaptic signaling currents (mEPSCs) in 
the nucleus accumbens shell and found a de- 
crease in frequency, but not amplitude, of mEPSCs 
(Fig. 4F). Lowered frequency of mEPSCs is a 
plausible functional consequence of reduced den- 
dritic spines, although we cannot rule out other 
adaptations. 

Our work has three major imphcations. First, 
we identify a mutation at the nucleotide level 
by using QTL analysis and provide evidence that 
CYFIP2 is a regulator of acute and sensitized co- 
caine response. Second, with over 20 commer- 
cially available C57BL/6 substrains, we describe 
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Fig. 4. CYFIP2 S968F causes destabilization, and knockout analyses of 
CyfipZ show that heterozygous mice have rescued cocaine response 
phenotypes. (A) Protein stability assay shows CYFIPZ S968F is less stable than 
wild type (WT). Human 293T cells transfected with FLAG-CYFIP2 (green bands) were 
treated cycloheximide for the indicated times, and Western blot was performed 
(bottom), p-tubulin was used as a control (red bands). Replicate experiments were 
quantitated, and half-life was calculated by using nonlinear one phase decay. K, the 
half-life parameter, is highly significant P < 0.0001 (top). (B) Quantitative IP 
experiments show CYFIPZ S968F interacts with WAVE complex members. IP (left) 
was conducted in replicate, and binding was quantitated (right). (C) Mice generated 
from ES cells harboring the knockout (KO) first allele of CyfipZ from the I KMC are 
on B6N background (designated Cyfipf^^^^'^ for WT and Cyppf^^- for 
heterozygous knockout). The homozygous deletion is lethal at birth and therefore 
not shown. Cyfip2^'^'^ animals have lower cocaine response than the Cyfip2^^'^~ 



mice. (D) The Cyfip2 heterozygous mice show higher sensitized response to cocaine 
than WT. Cyfipf^^^^^ (orange) and Cyfipf^'^- (purple) were injected with saline 
(S) or cocaine (C, 10 mg/kg). Painwise comparisons for genotype for each day were 
significant as indicated. (E) Diolistic labeling of nucleus accumbens neurons shows 
decrease in dendritic spine density in B6N. Sample images are shown with 
representative location of neurons (yellow on coronal section of brain). 
Quantitation and classification of spines are shown. At least six animals were used 
in each group with over 10 images collected from each animal. (F) B6N has a 
decrease in frequency of a-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid 
(AMPA) receptor mEPSCs. Sample traces of mEPSCs in nucleus accumbens shell 
medium spiny neurons from 86] and B6N mice. Dot plots for mEPSC amplitude 
and frequency. Data are mean ± SEM and Student's t test (*P < 0.05, **P < 
0.01, < 0.001, < 0.0001). Two-way ANOVA with Tukey post hoc test 
was used for (A). 



a framework to use mouse substrains as a rich 
genetic source for identifying previously unknown 
genes and alleles that regulate behavior. If other 
substrains such as C3H/He and DBA/2 are included, 
there may be hundreds of mouse substrains avail- 
able for analysis. Third, care must be taken when 
comparing new behavioral data from C57BL/6N 
substrains with existing data from C57BL/6J. 
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